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I.  INTRODUCTION 


The  JANNAF  Propulsion  Systems  Hazards  Subcommittee  s[)onsored  a 
Workshop  entitled  "Model  Predictions  of  the  Piston-Driven-Compaction  Experiment" 
'.vhich  was  held  on  25-26  October  1988  in  conjunction  with  the  25th  JANNAF 
Combustion  Meeting  at  NASA/Marshall  Space  Flight  Center,  Huntsville,  AL.  The 
general  problem  of  interest  is  the  transition  to  detonation  in  confined  granular 
energetic  material.  This  workshop  is  the  third  in  a  series  of  JANNAF  PSHS 
workshops  which  have  attempted  to  examine  various  aspects  of  the  problem.  The 
first  one  w'as  held  on  11-12  July  1985  in  the  Albuquerque  Convention  Center, 
Albuquerque  NM,  immediately  preceding  the  Eighth  Symposium  (International)  on 
Detonation,  15-19  July  1985.  The  objective  was  to  assess  our  knowledge  (or 
ignorance)  of  the  mechanical  properties  exhibited  by  compacted  granular  material. 
Discussion  concentrated  on  quasi-static  compaction,  as  well  as  dynamic  compaction 
and  small  amplitude  wave  propagation  which  may  involve  strain  rate  effects. 
Reference  [i]  is  the  workshop  summary. 

The  second  workshop  was  held  on  23  October  1986  at  the  NASA  /  Langley 
Research  Center  in  Hampton  VA,  in  conjunction  with  the  23rd  JANNAF 
Combustion  Meeting,  20-24  October  1986.  The  specific  objective  was  to  gather 
eight  modelers  together  and  compare  predictions  to  three  pre-assigned  "homework" 
problems  which  were  based  on  constant-velocity-piston  compaction  of  a  column  of 
confined  granular  material.  The  broader  objective  was  an  elementary  assessment 
of  the  predictive  capability  of  a  number  of  computer  codes  used  in  the  modeling 
of  Deflagration-to-Detonation  Transition  (DDT).  Of  course,  success  in  predicting 
compaction  wave  speed  is  only  a  neccessary  condition  and  not  sufficient  proof 
that  a  DDT  prediction  will  be  accurate.  The  good  news  from  this  workshop  was 
that  most  of  the  predictions  were  close  to  each  other  and  to  the  "correct" 
answers  which  were  taken  from  experimental  compaction  wave  speed  data 
obtained  by  H.  W.  Sandusky  (NSWC/WO)  in  his  Piston-Driven-Compaction  (PDC) 
apparatus.  The  bad  news  was  that  compaction  wave  speed  is  a  rather  forgiving 
parameter;  models  which  predict  nearly  the  same  wave  speed  were  observed  to 
predict  different  values  of  pressure  (or  stress)  in  the  aggregate.  It  was  concluded 
that  on  the  basis  of  compaction  wave  speed,  there  is  no  way  to  distinguish  the 
influence  of  (a)  strain-rate  sensitivity  during  dynamic  compaction,  (b)  compaction 
wave-induced  reaction  of  solid  material,  and  (c)  friction  between  the  sliding 

aggregate  and  the  wall  boundary.  Reference  [2]  is  the  workshop  summary. 

The  current  workshop  (Oct  1988)  is  an  attempt  to  confront  a  problem  closer 
to  "transition",  without  requiring  computations  of  a  full  high  order  detonation. 

The  aim  is  to  focus  on  wave  behavior  in  the  "intermediate  regime",  as  the  strong 
compaction  wave  accelerates  to  speeds  in  the  range  of  low  velocity  detonation. 
When  constructing  a  modeling  exercise  which  is  to  predict  the  behavior  of  an 
experiment,  an  absolute  necessity  is  the  ability  to  uniquely  specify  the  boundary 

conditions  and/or  source  terms  which  drive  the  problem.  If  the  problem  is  posed 
such  that  various  combinations  of  boundary  values  and  source  terms  could  be 
selected  as  reasonable,  then  comparison  of  model  predictions  will  be  meaningless. 
This  is  the  difficulty  associated  with  experiments  based  on  convective  ignitioii  - 

specifying  the  spatial  distribution  and  time-history  of  mass  flux  and  energy 
content  of  the  ignitor  system.  It  was  felt  that  Sandusky's  PDC  experiment 

should  not  suffer  from  tms  ambiguity;  the  experimental  values  of  projectile 

location  as  a  function  of  time  would  provide  a  unique  time-dependent  boundary 

condition  for  the  models.  Different  initial  wave  strengths  could  be  generated  by 
varying  the  initial  speed  of  the  projectile.  However,  the  original  PDC  experiment 
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would  have  to  be  modified  by  the  addition  of  one  or  more  pressure  tranoducers 
at  fixed  locations  along  the  tube  length.  These  would  provide  time-history  of 
stress  in  the  compacted  aggregate  as  the  event  developed,  and  help  distinguish 
among  various  model  predictions.  Thus,  the  current  workshop  extends  the  theme 
of  the  Oct  1986  workshop  [2]  to  a  more  difficult  problem. 

A  novel  aspect  of  the  current  workshop  is  that  the  modelers  were  asked  to 
make  true  predictions,  not  just  simulate  known  data.  Although  the  PDC 
experiments  were  run  in  advance,  the  results  were  not  released  until  the  second 
day  of  the  workshop.  Thus  all  the  modelers  (including  the  chairman)  were 
"blind".  Approximately  four  months  in  advance,  each  modeler  was  given  the 
initial  conditions,  the  time-dependent  boundary  condition,  and  properties  of  the 
granular  propellant.  There  were  four  experiments  involving  two  different 

materials.  The  modeler  was  asked  to  predict  the  compaction  wave  locus,  and  the 
stress  time-history  which  would  be  seen  by  gages  at  two  fixed  locations  in  the 
chamber.  Six  individuals  agreed  to  be  subjected  to  this  torture,  and  five  brought 
predictions  to  the  Workshop.  They  were: 

M.  R.  Baer,  Sandia  National  Laboratory  /  Albuquerque 

K.  Kim  and  T.  Hsieh,  Naval  Surface  Warfare  Center  /  White  Oak 

Laboratory 

D.  E.  Kooker,  U.  S.  Army  Ballistic  Research  Laboratory 

C.  F.  Price,  Naval  Weapons  Center  /  China  Lake 

A.  M.  Weston,  Brobeck  Corp  /  Lawrence  Livermore  National  Laboratory 

The  first  day  the  workshop  was  devoted  to  presentation  of  the  models,  and 
comparison  of  the  various  predictions  with  each  other.  H.  W.  Sandusky 

(NSWC/WO)  unveiled  the  experimental  data  on  the  second  day,  and  the  remainder 
of  the  time  was  spent  comparing  model  predictions  with  this  data. 


[Experiment  by  Sandusky  et.  al./NSWC] 

Wall  Pressure  (Stress)  -  Time  History 


76.2  mm 


38.1  mm 


'Mum  or  Steel  Tube 


- s 


Granular  Propellant 


Given; 


Ionization  Pins  or  Self-Shorting  Probes 

(a)  propellant  properties  - 

QS  compaction,  grain  size,  burning  rate,  etc. 

(b)  projectile  time-history 

Predict  Time-History:  (I)  Compaction  Front  Location 

(2)  Pressure  at  Two  Fixed  Locations 


Fig.  1  -  Schematic  of  Modeling  Problem  and  NSWC  Experiment  [3]. 
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II.  PROBLEM  DESCRIPTION 


The  modeling  exercise,  shown  schematically  in  Fig.  1,  is  to  simulate 
Sandusky's  Piston-Driven-Compaction  (PDC)  experiment  in  which  a  quiescent 
granular  sample  is  confined  in  an  aluminum  or  heavy-wall  steel  tube  (25.4  mm 

ID),  and  then  subjected  to  impact  by  a  cylindrical  blunt-face  Lexan  projectile 
(300  mm  long  by  25.4  mm  in  diameter).  Diagnostics  include: 

(a)  ionization  pins  and/or  self -shorting  probes  at  uniform  intervals  along  the 
length  of  the  tube, 

(b)  microwave  interferometry  through  the  downstream  end  of  the  tube 

(looking  toward  the  impacted  end),  and 

(c)  flush-mounted  pressure  transducers  in  the  tube  wall  at  38.1  mm  and  76.2 
mm  from  the  impacted  end. 

The  combination  of  (a)  and  (b)  serves  to  locate  the  leading  compaction  front  as  a 
function  of  time.  Further,  the  location  of  the  projectile  face  after  impact  can  be 
determined  by  the  position  of  scribe  lines  on  the  cylindrical  projectile 

photographed  through  a  slit  in  the  tube  wall  near  the  impacted  end;  this  provides 
the  time-dependent  boundary  condition  for  the  modelers.  More  detail  about  this 
series  of  experiments  can  be  found  in  Sandusky  et.  al.  [3].  As  the  basis  for  this 
workshop,  there  were  four  PDC  runs  involving  two  different  ball  propellants.  Two 
primary  reasons  motivated  the  selection  of  these  materials;  (1)  compaction  of  an 
aggregate  involves  plastic  deformation  of  grains  [4],  but  not  fracture  [5],  and  (2) 
they  are  nearly  spherical  particles  which  might  be  expected  to  form  a  regular 

lattice  which  may  be  amenable  to  analysis  [6].  Properties  of  both  propellants  are 
listed  in  Table  I. 

TABLE  I.  -  PROPELLANT  PROPERTIES 


Propellant  Name 

WC-140 

TS-3659 

Composition: 

NC  (nitrogen  content  =  13.05%) 

98.53% 

76.98% 

NG 

0.00% 

21.57% 

Stabi 1 izer/Additive 

1.47% 

-- 

Diphenylamine 

— 

1.14% 

Dinitrotoluene 

— 

0.31% 

Theoretical  Maximum  Density  (g/cc) 

1.65 

1.64 

Heat  of  Explosion  (cal/g,  experimental) 

919 

1104 

Nominal  Particle  Size  (/:zm) 

411 

434 

(spheres  &  "jelly  beans") 

(spheres) 
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Information  Supplied  to  Modelers 

(A)  From  the  PDC  Experiment: 

The  initial  conditions  for  the  modeling  exercises  were  taken  directly  from  the 
N'SWC  experiments  [3],  Listed  in  Table  II  below  are  values  for  tube  length  (L), 
initial  density  (9fcTMD^),  the  projectile  velocity  (V  )  just  before  impact,  and  the 

initial  particle  velocity  (u)  and  compaction  wave  speed  (U)  just  after  impact.  Of 
course,  the  values  of  particle  velocity  (the  velocity  of  the  projectile  face  after 
impact)  and  compaction  wave  speed  will  not  necessarily  remain  constant. 
However,  the  modeler  was  also  given  the  location  of  the  projectile  as  a  function 
of  time  after  impact  which  could  be  used  as  a  "piston"  boundary  condition;  this 
avoids  all  the  uncertainties  associated  with  a  description  of  the  actual  impact 
event. 


Sandusky's  [7^  original  PDC  experiment  performed  in  a  Lexan  tube  produced 
the  correlation  "t  At  =  constant",  where  t  is  a  measure  of  the  stress  in  the 
compacted  solid  phase  and  At  is  the  time  interval  after  impact  before  the 
appearance  of  visible  flame  (usually  in  material  adjacent  to  the  projectile  face). 
It  was  intended  that  this  modeling  exercise  w'ould  provide  values  of  At  for  each 
propellant,  which  the  the  modelers  could  use  to  help  calibrate  parameters  in  their 
reaction  models.  Several  PDC  runs  were  made  in  a  Lexan  tube  w'ith  TS-.3659 
propellant,  but  the  camera  records  showed  no  evidence  of  visible  light;  other 
observations,  however,  clearly  indicated  that  reaction  was  underwav  This 
"negative"  finding  may  have  important  implications  concerning  the  extent  and/or 
type  of  reaction  which  begins  the  event.  As  a  substitute  in  the  four  modeling 
problems,  Sandusky  examined  both  gage  records  and  estimated  the  time  at  which 
rapid  pressurization  occurred.  These  two  points  in  the  x-t  plane  were  used  to 
linearly  extrapolate  back  to  the  location  of  the  projectile  face.  The  value  of  time 
at  the  intersection  defined  a  "At",  the  time  interval  before  the  onset  of  rapid 
pressurization  (not  the  appearance  of  visible  flame).  The  modelers  were  provided 
a  value  (in  one  case,  a  range  of  values)  of  "At"  for  each  of  the  problems. 

TABLE  11.  -  INITIAL  CONDITIONS  FOR  PDC  MODELING  PROBLEMS 


PROPELLANT:  TS-3659 


Problem 

(PDC  name) 

%TMDq 

L(mm) 

Tube 

Vp(m/s) 

u(m/s) 

U(m/s) 

#1 

PDC-80 

60.2 

101.7 

Alum 

160 

127 

494 

#2 

PDC-81 

60.1 

146.8 

Steel 

237 

192 

534 

#3 

PDC-82 

60.1 

146.8 

Steel 

300 

216 

557 

PROPELLANT:  WC-140 

#4 

PnC-M30 

60.5 

146.8 

Steel 

ZIO 

161 

563 

^  .t  4c  :4c  4c  34c  4c  4=  4^  4c  *  4c  ^  4: ^  ^  4c  4c  ]|c  4c  %  4c  4c  4c  ^  4;  4c  4c  :4c  4:  ie  sfc  ^  ^  4c  :4c  *  4c  t  :4c  4c  :4c  :|c  4c  %  4: 

♦  ♦*ir****+4c4c*4c**4c*******4c*4*4*****+4*****4*4*****4t*****4*4*4444c*44c444c444c444+444c4c4c4c4c+**4c4c**44c4c+*4c 
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(B)  Quasi-Static  Compaction  Data; 

Current  models  of  compacted  granular  material  are  not  yet  capable  of 
predicting  the  state  of  stress  (including  porosity)  in  the  compressed  aggregate.  All 
of  the  theories  involve  parameters  which  must  be  calibrated  by  matching 

experimental  data.  For  this  purpose,  quasi-static  compaction  experiments  were 
performed  over  a  range  of  mixture  densities  on  both  ball  propellants  using  the 
double-piston  in  cylinder  apparatus  discussed  in  Campbell,  Elban  and  Coyne  [8]. 
Data  supplied  to  the  modelers  included  axial  stress,  radial  stress,  force  applied  to 
the  driving  ram,  force  transmitted  to  the  passive  ram,  and  the  length  of  the 
sample  volume.  Differences  between  applied  and  transmitted  force  (along  with  the 
length  measurement)  allowed  computation  of  the  wall  friction  force.  Determination 
of  porosity  required  the  modeler  to  estimate  the  compressibility  of  the  solid 

material.  For  more  detail  on  the  experiment,  see  Sandusky  et.  al.  (3). 

(C)  Propellant  Ignition  and  Burning  Rate  Data: 

A  reasonable  description  of  the  burning  rate  of  single-base  WC-140 
propellant  over  the  pressure  range  100  psi  -  25  Kpsi  was  provided  by  a 
combination  of  recent  strand  burner  data  obtained  by  Price  and  Atwood  [9] 
supplemented  with  BRL  [10]  closed-bomb  data.  Reference  [9]  also  made  available 
some  recent  ignition  data  generated  from  both  arc  image  and  CO2  laser  sources. 
For  the  double-base  propellant  TS-3659,  burning  rate  information  was  incomplete. 
Strand  burner  data  were  nonexistent.  This  exercise  relied  entirely  on  NWC/China 
Lake  closed-bomb  data  [11].  However,  as  is  typical  of  the  closed  bomb,  rates 
below  approximately  2-3  Kpsi  show  the  anomolous  "tail  off".  Without  low 
pressure  strand  data,  the  modelers  were  left  to  guess  the  rates  at  lower  pressure, 
using  the  behavior  of  WC-140  as  a  guide.  At  the  moment,  there  are  no 
estimates  of  how  sensitive  the  wave  acceleration  process  is  to  variations  in  low 
pressure  burning  rate. 

(D)  Compressibility  of  TMD  Propellant: 

Although  Hugoniot  curves  have  been  determined  for  many  inert  and  explosive 
compounds,  no  information  could  be  found  for  either  WC-140  or  TS-3659 
propellant.  As  a  substitute,  some  Russian  data  [12]  for  N-powder  with  a 
composition  of  approximately  40%  NG  in  NC  was  provided,  viz., 

V  (km/s)  =  1.7  4-  1.85  *  u  (km/s). 
s  s 

With  40%  NG,  this  material  will  be  slightly  softer  than  either  WC-140  or 
TS-3659.  but  the  influence  on  the  present  modeling  computations  should  be 
minimal. 

Given  this  information,  the  modeler  was  asked  in  each  case  to  provide 
predictions  for  (1)  compaction  wave  locus,  and  (2)  the  time-history  of  whatever 
stress  component  a  wall-mounted  transducer  should  respond  to  at  both  gage 
locations.  Because  of  complications  involved  with  wave  reflection  from  an  end 
plug  which  was  not  well  characterized,  it  was  agreed  that  the  computational 
problem  will  terminate  when  the  leading  compaction  front  strikes  the  downstream 
boundary  (the  far  end  of  the  PDC  tube). 
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III.  DISCUSSION  OF  RESULTS 


Although  the  five  models  employed  here  were  designed  to  solve  a  similar 
problem,  there  are  many  important  differences  which  lead  to  a  diversity  of 
predictions.  Unfortunately,  the  space  limitations  of  this  workshop  summary  will 
not  permit  a  detailed  discussion  of  the  assumptions,  mathematical  formulation,  and 
numerical  solution  techniques  in  each  model.  To  help  orient  the  reader,  the 
Appendix  contains  a  schematic  overview  of  important  aspects  of  each  model. 
These  schematics  are  not  intended  as  a  complete  picture  of  the  model;  the  cited 
references  should  be  consulted  for  further  detail,  or  questions  should  be  addressed 
to  the  individual  modeler  in  cases  where  documentation  is  not  yet  available.  In 
the  figures  below  which  show  comparison  of  model  predictions  to  the  experimental 
data  of  Sandusky  et.  al.  [3],  not  all  of  the  models  are  represented  on  each  figure. 
This  occurs  because  not  all  of  the  modelers  brought  predictions  for  all  four 
exercises.  Table  III  below  can  be  used  as  a  guide. 

TABLE  III.  -  SUMMARY  OF  CONTRIBUTIONS 

JANNAF  PSHS  WORKSHOP 
"Model  Predictions  of  the  PDC  Experiment" 

25-26  October  1988,  NASA/MSFC,  Huntsville,  AL 


PREDICTIONS  FOR- 
CONTRIBUTOR  | 

PDC-80 

#1 

PDC-81 

#2 

PDC-82 

#3 

PDC-M30 

Baer  /  SNLA 

Yes 

Yes 

Yes 

Yes 

Kim  Sc  Hsieh/NSWC 

Yes 

Yes 

Kooker  /  BRL 

Yes 

Yes 

Yes 

Yes 

Price  /  NWC 

Yes 

Yes 

Yes 

Yes 

Weston  /  LLNL 

Yes 

*********************************************************************************************** 

Figures  2  through  9  show  the  comparison  of  modeling  predictions  with  the 
experimental  data.  The  figures  are  arranged  in  pairs;  Fig.  2  is  the  compaction 
wave  locus  and  Fig.  3  contains  the  two  transducer  time-histories  for  PDC-80 
(Problem  #1).  Figures  4  and  5  pertain  to  PDC-81  (Problem  #2),  etc.  In  each 
figure,  the  top  portion  (a)  is  only  the  experimental  data,  and  the  bottom  portion 
(b)  shows  the  comparison  of  the  model  predictions  with  that  data.  Many  of  the 
comparison  figures  ((b)  portion]  are  busy  and  somewhat  confusing.  A  vertical  line 
with  the  notation  "end  of  problem"  appears  on  each  figure,  and  denotes  the  value 
of  time  at  which  the  experimental  data  indicate  the  leading  compaction  front 
arrived  at  the  end  of  the  PDC  tube.  Beyond  this  value  of  time,  the 
computational  problem  is  over. 
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The  nature  of  this  problem  is  deceptive.  The  initial  impact  generates  a 

strong  compaction  wave  which  propagates  into  the  quiescent  bed.  Taking  a 

simplistic  view,  not  niuch  happens  behind  this  wave  for  a  certain  time  interval; 

then  suddenly  the  pressure  can  rise  dramatically.  The  deception  is  in  the  region 
behind  the  compaction  front  where  the  pressure  transducers  may  indicate  that  not 
much  is  happening.  In  fact  there  is  apparently  an  intense  competition  between 
sources  of  heat  generation  (reactions)  and  heat  loss  mechanisms.  The  outcome  of 
the  competition  produces  the  "timing"  of  a  runaway  (if  one  occurs).  On  the  basis 
of  the  modeling  predictions,  a  small  change  in  a  reaction  parameter  can  lead  to  a 
major  change  in  the  timing  of  the  runaway.  Probably  as  no  surprise,  none  of 
the  models  got  everything  "correct".  A  major  conclusion  from  the  Oct  1986 
Workshop  [2]  is  reinforced  by  the  results  from  the  present  comparisons; 

compaction  wave  speed  (and  hence  location)  is  a  forgiving  parameter.  A  close 
match  of  compaction  wave  locus  does  not  necessarily  mean  a  close  match  with 
the  pressure  transducers, 

PDC-80  (Pj;ob]em  #1) 

PDC-80,  based  on  TS-3659  propellant,  was  run  in  the  short  tube  (101.7  mm) 
and  had  the  slowest  initial  projectile  speed  (160  m/s).  As  shown  in  Fig.  2a,  the 
compaction  wave  locus  is  determined  by  two  points  dictated  by  the  reporting  time 
from  each  pressure  transducer;  PDC-80  was  the  only  run  which  did  not  have 
microwave  interferometry.  As  shown  in  Fig.  2b,  alf  three  predictions  are  close  to 
the  first  data  point,  with  some  divergence  toward  the  latter  stages  of  propagation. 
The  sudden  "turn  up"  predicted  by  Baer  (SNLA)  cannot  be  confirmed  with  the 
two  data  points  available.  The  "slow  down"  predicted  by  Price  (NWC)  may  be 
partly  the  result  of  his  time-dependent  boundary  condition.  In  all  his 
computations.  Price  (NWC)  elected  to  use  a  lumped-mass  solution  of  projectile 
impact  on  the  granular  bed,  rather  than  employ  the  given  experimental  projectile 
location  as  a  function  of  time  as  a  "piston"  boundary  condition.  This  solution 
predicted  a  time-dependent  location  of  the  projectile  face  which  differed  from  the 
experimental  observation. 

The  time-history  from  both  pressure  transducers  is  shown  in  Fig.  3a.  Both 
gages  confirm  the  onset  of  rapid  pressurization,  although  the  rapid  rise  of  gage 
#2  occurs  beyond  the  time  (~204  tts)  at  which  the  leading  compaction  front 
reached  the  end  of  the  PDC  tube.  There  are  two  common  features  in  the  gage 

responses  from  all  four  PDC  runs.  The  first  one  is  clearly  evident  here  in  the 

record  from  Gage  #2,  and  to  a  lesser  extent  from  Gage  #1.  This  is  the  almost 
linear  "decay"  of  the  report  value  after  the  gage  has  been  uncovered,  before  the 
beginning  of  rapid  pressurization  presumeably  due  to  significant  reaction  [see  Fig. 
10  also].  None  of  ^  models  predicted  this  behavior.  There  was  considerable 

discussion  about  the  possible  physics  involved,  as  well  as  the  possibility  of  gage 
anomalies.  One  suggestion  was  rate-dependent  stress  behavior  in  the  suddenly 
compacted  aggregate,  which  might  influence  the  deviator  stress  component  and 
hence  play  a  role  in  the  stiess  value  which  is  being  measured  by  the 
wall-mounted  transducer.  The  issue  was  not  resolved.  The  second  feature  is  the 
smaller  magnitude  of  initial  response  (and  subsequent  values  before  reaction 
begins)  reported  by  Gage  #2  compared  to  Gage  #1.  This  is  most  likely  the 
influence  of  wall  boundary  friction.  Of  the  three  model  calculations  illustrated  in 
Fig.  3b,  Price's  (NWC)  model  was  the  only  one  to  account  for  wall  friction  and  it 
predicts  this  effect.  Figure  3a  also  illustrates  another  trend  when  comparing 
model  predictions  with  the  transducer  records.  Both  Kooker  (BRL)  and  Baer 

(SNLA)  predict  an  initial  level  of  gage  response  which  exceeds  the  experimental 
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t  (fJS  ) 


Fig.  2a  -  Compaction  Wave  Locus  for  PDC-80  (ProWem  #1);  160  m/s  Impact  on  60.296  TMD  TS-3659. 
Tube  Length  is  101.7  mm.  Experimental  Data  from  Sandusky  et  al.  [3). 


Fig.  2b  -  Compaction  Wave  Locus  for  PDC-80  (Problem  #1);  160  m/s  Impact  on  60.296  TMD  TS-3659. 
Comparison  of  Workshop  Predictions  with  Data  from  Sandusky  et.  al.  (3|. 
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Fig.  3a  -  Time-History  From  Wall-Mounted  Transducers  for  PDC-80  (Problem  #1);  160  m/s  Impact  on 
60.296  TMD  TS-3659.  Experimental  Data  from  Sandusky  et  al.  (3).  Gage  #1  is  Located  at  38.1  mm, 
and  Gage  #2  at  76.2  mm;  Tube  Length  is  101.7  mm. 


Fig.  3b  -  Time-History  From  Wall-Mounted  Transducers  for  PDC-80  (I’robletii  #1);  160  m/s  Impact  on 
60.296  TMD  TS-3659.  Comparison  of  Workshop  Predictions  with  Data  from  Sandusky  et  al.  |3|.  Gage 
#1  is  Located  at  38.1  mm,  and  Gage  #2  at  762  mm;  Tube  Length  is  101.7  mm. 
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value.  The  mathematical  formulation  in  all  three  models  (at  the  time  of  the 
Workshop)  treated  the  stress  tensor  of  the  compressed  aggregate  as  a  pressure, 
and  ignored  deviator  components.  Both  Kooker  (BRL)  and  Baer  (SNLA)  based 
their  predictions  for  gage  response  on  this  value  of  mixture  pressure.  However, 
the  quasi-static  compaction  data  (S)  for  TS-3659  distinctly  shows  the  presence  of 

a  deviator  component  which  leads  to  a  significant  reduction  in  the  value  of  radial 
stress.  The  data  [3]  indicate  that  the  ratio  of  the  radial  stress  component  to  the 

axial  component  varies  with  the  state  of  compaction.  Price  (NWC)  chose  an 

average  value  of  0.6  as  representative  of  the  data,  and  multiplied  his  predictions 

for  pressure  by  this  value  to  get  gage  response.  As  the  comparison  in  Fig.  3b 

indicates,  these  values  are  more  realistic.  The  important  implication  is  that  the 
dynamically-formed  aggregate  in  both  materials  apparently  supports  deviator  s..ress 
components  which  are  similar  to  the  quasi-static  stress  state.  The  wall-mounted 
transducers  are  responding  to  the  radial  component. 

Compaction  wave  thickness  can  also  be  estimated  from  the  transducer 

records,  assuming  wave  speed  is  known.  For  example,  Gage  #1  in  Fig.  3 

indicates  a  rise  time  of  approximately  12  iis  which  implies  a  wave  thickness  of  6 
mm  if  the  speed  is  taken  as  0.5  mm/Ms.  Recall  that  Sandusky's  [7]  earlier  work 
utilizing  flash  radiography  in  a  Lexan  tube  estimated  the  wave  thickness  as  the 

distance  between  two  tracer  wires,  i.e.,  6.35  mm.  Model  predictions  of  wave 
thickness  (or  gage  rise  time)  are  significantly  influenced  by  the  equation 
formulation  and  the  numerical  solution  technique.  If  the  equation  system  neglects 
the  diffusive  effects  of  Fourier  heat  conduction  (not  heat  transfer  between  phases), 
viscous  stresses  (not  the  drag  force  between  phases),  and  rate-dependent  resistance 
forces  during  porosity  adjustments  in  the  aggregate  (often  modeled  with  a 
"compaction  viscosity"),  then  the  mathematical  solution  for  a  compaction  wave  is  a 
shock  wave  with  zero  thickness.  If  the  numerical  solucion  predicts  otherwise,  the 
finite  thickness  is  due  to  artificial  diffusion  added  by  the  numerical  integration 
technique.  Baer's  (SNLA)  formulation,  which  does  account  for  rate-dependent 
porosity  adjustments  and  hence  includes  a  physical  basis  for  a  finite  thickness 
wave,  predicts  a  gage  rise  behavior  quite  similar  to  the  experimental  data.  This 
is  a  definite  achievement.  Kooker  (BRL)  assumes  instantaneous  adjustment  (zero 
compaction  viscosity)  to  the  equilibrium  stress  state  and  then  tracks  compaction 
waves  as  infinitesimally-thin  shock  waves  with  the  method-of-characteristics;  the 
predictions  show  a  zero-thickness  wave  which  passes  the  gage  location  near  the 
mid-point  of  the  rise  time.  Price’s  (NWC)  formulation  also  neglects  rate- 

dependent  resistance  (other  than  inertial  forces)  to  porosity  change,  but  the 
solution  predicts  a  wave  thickness  somewhat  greater  than  the  data. 

Both  gage  records  in  Fig.  3  confirm  that  some  type  of  reaction  process  was 
under  way  near  the  end  of  this  PDC  experiment.  The  comparisons  in  Fig.  3b 
also  testify  to  the  difficulty  the  models  had  in  predicting  the  timing  of  this 

pressure  rise.  The  input  data  supplied  to  the  modelers  was  insufficient  to 

calibrate  all  the  "reaction  parameters"  in  the  various  formulations.  The  remaining 
values  had  to  be  "guessed",  or  estimated  by  comparison  to  other  available  data  on 
similar  materials.  Unfortunately,  a  slight  variation  in  some  of  these  values  could 
often  mean  the  difference  between  benign  behavior  and  a  runaway.  As  seen  in 
Fig.  3b,  Baer's  (SNLA)  computations  show  a  premature  runaway,  Kooker's  (BRL) 

computations  show  a  turn-up  closer  to  the  gage  records  (probably  fortuitous),  and 
Price's  (NWC)  computations  indicate  no  runaway. 
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PpC-81  (Problem  #2) 

PDC-81  is  also  based  on  TS-3659  propellant,  but  was  run  in  the  146.8  mm 
long  heavy-wall  steel  tube  using  a  higher  initial  projectile  speed  of  237  m/s.  The 
leading  compaction  front  locus,  as  shown  in  Fig.  4a,  is  defined  in  detail  by  the 
data  from  the  microwave  interferometer.  Compaction  wave  speed  is  virtually 
constant  at  534  m/s  until  late  in  the  run  when  it  abruptly  increases  to  1800  m/s. 
Comparing  model  predictions  to  this  data  in  Fig.  4b  shows  that  all  four 
calculations  are  close  to  the  experiment  until  approximately  half  way  down  the 
tube.  Baer's  (SNLA)  model  prematurely  predicts  an  abrupt  increase  in  speed,  but 
the  new  velocity  is  nearly  equal  to  what  the  experiment  will  eventually  do.  This 
is  encouraging.  Again,  timing  of  the  event  can  be  difficult  to  predict.  Kim  and 
Hsieh's  (NSWC)  model  predicts  a  violent  runaway  event  with  the  wave  about 
midway  down  the  tube,  and  this  halts  the  calculation  (see  the  predicted  pressure 
history  in  Fig.  5b).  Kooker's  (BRL)  model  predicts  a  mild  turn-up  just  as  the 
wave  approaches  the  end  of  the  tube,  and  Price's  (NWC)  model  does  not  yet 

show  the  effects  of  reaction. 

The  model  predictions  for  wall  stress  time-history  are  compared  to  the 

transducer  records  in  Fig.  5.  For  each  of  the  four  models,  the  predicted  behavior 
of  the  compaction  front  locus  discussed  above  is  mirrored  in  their  stress 
predictions.  Although  the  overall  slopes  of  a  couple  of  the  runaway  curves  are 
similar  to  the  data,  the  predicted  timing  of  the  event  is  poor.  Some  discussion 
in  the  workshop  arose  over  the  "plateau"  behavior  (near  200  MPa)  of  the  Gage 
#1  response  between  180-210  (is.  The  question  was  what  mechanism  is  resposible 
for  the  20-30  (is  "pause"  in  the  growth  of  pressure?  One  speculation  was  an 

anomaly  in  "hot  spot"  compressive  combustion.  Another  speculation  was  the 
influence  of  expansion  waves  generated  by  a  slow  down  (or  reversal  of  direction?) 
of  the  projectile,  i.e.,  a  loss  of  confinement.  This  second  speculation  highlights  a 
nagging  problem  in  these  exercises.  Because  of  experimental  difficulties,  the 

position  of  the  projectile,  after  impact,  is  sometimes  known  only  during  the  early 
time  of  the  experiment  (in  PDC-81,  only  up  to  92  (is).  Furthermore,  in  those 
runs  based  on  higher  projectile  impact  speeds,  there  is  little  doubt  that  the  Lexan 
projectile  deforms;  Sandusky  attempts  to  account  for  this  deformation  when 
interpreting  the  position  of  the  scribe  lines  on  the  body  of  the  projectile.  Thus, 
(a)  the  modeler  may  not  have  the  complete  time-dependent  boundary  condition, 
and  (b)  its  accuracy  is  in  question.  Sandusky  suggested  that  any  future 
experiments  of  this  type  should  be  run  with  an  aluminum  projectile. 

P^-82  (Problem  #3) 

PDC-82  is  still  based  on  TS-3659  propellant  confined  in  the  146.8  mm  long 
heavy-wall  steel  tube,  but  involved  a  300  m/s  impact.  This  was  a  more 

energetic  event  which  produced  an  increase  in  speed  to  2110  m/s  as  shown  in  Fig. 

6  by  the  compaction  front  locus  defined  by  the  microwave  interferometry.  The 

model  predictions  for  compaction  front  locus  are  fairly  close  to  the  data  until  the 
vicinity  of  the  turn.  The  two  closest  predictions  were  Price's  (NWC)  model  which 
turned  up  early  and  Baer's  (SNLA)  model  which  was  slightly  late,  but  both 

suggested  a  speed  in  the  correct  range.  Kim  and  Hsieh's  (NSWC)  model  did  not 
show  a  significant  change  in  wave  speed,  but  the  pressure  field  is  increasing 
dramatically  at  the  end  of  the  computation.  Kooker's  (BRL)  model  ran  into 
numerical  difficulties  in  attempting  the  turn,  which  terminated  the  calculation. 

(text  continues  on  page  26). 
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Fig.  4a  -  Compaction  Wave  Locus  for  PDC-81  (Problem  #2);  237  n»/s  Impact  on  60.1%  TMD  TS-3659. 
Tube  Length  is  146.8  mm.  Experimental  Date  of  Sandusky  et,  al.  (3). 


Fig.  4b  -  Compaction  Wave  Locus  for  PDC-81  (Problem  #2);  237  m/s  Impact  on  60.196  TMD  TS-3659. 
Comparison  of  Workshop  Predictions  with  Microwave  Date  from  Sandusky  et.  al.  (3). 
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t  Ou) 

Fig.  5a  -  Time-History  From  Wall-Mounted  Transducers  for  PDC-81  (Problem  #2);  237  m/s  Impact  on 
60.19i  TMD  TS-3659.  Experimental  Data  from  Sandusky  et  aL  [3).  Gage  #1  is  Located  at  3U  mm, 
and  Gage  #2  at  76.2  mm;  Tube  Length  is  1463  mm. 


73 

I 


s 


—100  “I  I  I  I  I '  I  "  '  r  I  I  f  I  I  I 

0  25  50  75  100  125  150  175  200  225  250  275  300 

111116  /uLS 


Fig  5b  -  Time-History  From  Wall-Mounted  Transducers  for  PDC-81  (Problem  #2);  237  m/a  Impact  on 
60.19(  TMD  TS-3659.  Comparison  of  Workshop  Predictions  with  Data  from  Sandusky  et  al.  |3|.  Gage 
#1  is  Located  at  38.1  mm,  and  Gage  #2  at  763  mm;  Tube  Length  is  146.8  mm. 
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Fig.  6b  -  Compaction  Wave  Locus  for  PDC-82  (Problem  #.3);  .300  m/s  Impact  on  60.1%  TMU  TS-3659. 
Comparison  of  Workshop  Predictions  with  Microwave  Data  from  Sandusky  et.  al.  13]. 
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.  end  of  problem 
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Fig.  7a  -  Time-Mislory  From  Wall-Mounted  Transduc.ors  for  l’DC-82  (Problem  #3);  300  in/s  Impact 
on  60.1%  TMD  TS-3659.  Experimental  Data  from  Sandusky  et.  al.  (3].  Gage  #1  is  Located  at  38.1 
mm,  and  Gage  #2  at  76.2  mm;  Tube  Length  is  146.8  mm. 


Fig.  8a  -  Compaction  Wave  Locus  for  PDC-M30  (Problen.  #4);  210  m/s  Impact  on  60.5%  TMD 
WC-140.  Tube  Length  is  146.8  mm.  Experimental  Data  from  Sandusky  eb  al.  |3j. 
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Fig.  8b  -  Compaction  Wave  Locus  for  PDC-M30  (Problem  #4);  210  m/s  Impact  on  60.5% 
WC-140.  Comparison  of  Workshop  Predictions  with  Microwave  Data  from  Sandusky  et.  al.  |3|. 
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Fig.  9a  -  Time-History  From  Wall-Mounted  Transducers  for  PDC-M30  (Problem  #4);  210  m/s  Impact 
on  60.59(  TMD  WC-140.  Experimental  Data  from  Sandusky  et  al.  |3|.  Gage  #1  is  Locate  at  38.1 
mm,  and  Gage  #2  at  76.2  mm;  Tube  Length  is  146.8  mm. 


Tiine  /us 

Fig.  9b  -  Time-History  From  Wall-Mounted  Transducers  for  PDC-M30  (Problem  #4);  210  m/s  Impact 
on  60.59t  TMD  WC-MO.  Comparison  of  Workshop  Prrtliriinns  with  Dnin  from  Snndiisky  cL  nl.  |;i|. 
Gage  #1  is  Located  at  38.1  mm,  and  Gage  #2  at  70.2  mm;  Tube  Length  is  146J  mm. 
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The  transducer  records  from  PDC-82  shown  in  Fig.  7  confirm  the  strength 
of  the  combustion  event.  Gage  #2  sees  a  strong  shock-like  wave,  while  Gage  #1 
displays  some  unusual  undulations.  None  of  the  models  came  very  close  to  these 
records.  Possibly  projectile  motion  after  impact  has  led  to  a  complicated  system 
of  expansion  waves  which  influence  the  response  at  Gage  #1;  Sandusky  (NSWC) 
commented  that  the  reaction  in  this  run  was  violent  enough  to  virtually  force  the 
projectile  back  out  of  the  tube. 

PDC-xM30  (Problem  #4) 

PDC-M30  changed  propellant  to  the  less  energetic  single-base  WC-140,  but 
w'as  conducted  in  the  146.8  mm  long  heavy-wall  steel  tube.  Projectile  impact 

speed  was  210  m/s.  Compaction  front  locus  is  displayed  in  Fig.  8,  which  shows 
an  unusual  behavior  after  approximately  270  iis  as  the  compacted  bed  forces  the 
end  plug  out  of  the  tube.  The  microwave  interferometry  has  tracked  the  motion 
of  a  second  wave  front  (apparently  driven  by  combustion)  in  the  bed  as  it  is 
being  expelled.  As  a  conciliatory  gesture,  the  modelers  were  not  held  accountable 
for  this  late  time  behavior.  Up  to  270  tfs,  the  compaction  front  moves  at  nearly 
constant  speed,  which  most  of  the  models  predicted  until  the  onset  of  significant 
reaction.  Price  (NWC)  offered  two  predictions  which  differed  only  by  a  slight 
change  in  one  of  his  reaction  parameters;  the  influence  on  the  predicted  wave 
locus  was  substantial. 

The  transducer  time-histories  are  compared  to  the  model  predictions  in  Fig. 
9.  Unfortunately,  the  significant  rise  seen  by  both  gages  occurred  after  270  us 

when  the  compaction  front  reached  the  downstream  end  of  the  tube.  Baer 

(SNLA)  and  Kooker  (BRL)  were  close  to  the  reporting  time  of  the  gages,  while 
Price's  (NWC)  "number  1"  prediction  gives  a  better  representation  of  the  low 
magnitudes  from  both  transducers  {see  Fig.  10].  Weston's  (LLNL)  predictions  are 
compared  to  the  data  separately  in  Fig.  11  and  seem  to  show  evidence  of  an 
instability  (Weston  discovered  an  error  in  his  gas-phase  equation-of-state  after  the 

Workshop).  None  of  the  predictions  could  be  considered  as  accurate.  Again, 

timing  of  the  runaway  pressure  field  is  a  delicate  prediction. 

IV.  CONCLUDING  REMARKS 

During  the  course  of  the  Workshop,  a  number  of  discussions  arose  about  the 
behavior  of  the  strong  compaction  wave  and  how  the  theories  are  modeling  its 
role  in  initiating  reaction.  Apparently  there  are  many  cases  where  the 

experimental  value  of  compaction  wave  speed  exceeds  the  model  prediction  for  an 
inert  wave,  when  the  model  is  based  on  rate-independent  porosity  change  which  is 
calibrated  with  the  quasi-static  compaction  data.  The  question  of  what 

mechanism(s)  is  responsible  for  the  increase  in  speed  seems  to  have  provoked  two 
F^'hools  of  thought.  The  first  one  suggests  that  a  rate-dependent  visco-plastic 
material  response  of  the  aggregate  would  produce  a  stiffer  material  and  thus 
support  a  higher  wave  speed.  This  dynamic  stress  state  would  look  nearly 

time-invariant  to  the  PDC  experiment  if  the  aggregate's  characteristic  relaxation 
time  is  large  enough.  The  question  of  how  this  behavior  would  influence  the 

compaction  wave  thickness  has  not  been  ansv/ered.  The  question  of  how  the 

reaction  begins  in  the  compressed  aggregate  still  remains,  "rhe  second  school  of 
thought  is  that  these  compaction  fronts  are  not  inert;  the  high  strain  rate  in  the 
wave  front  itself  begins  the  reaction  process.  Of  course,  many  new  questions 
now  arise  concerning  what  type  of  reaction  begins  and  how  much  energy  is 

released.  There  is  little  consensus  from  this  point  on.  There  is  also  the 
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Fig.  10  -  Time-History  From  Wali-Mounted  Transducers  for  PDC-M30  (Problem  #4);  210  m/s  Impact 
on  60.5%  TMD  WC-140.  Comparison  of  Workshop  Predictions  with  Data  from  Sandusky  et.  al.  [3|. 
Scale  to  Illustrate  Low  Level  Gage  Response. 


t  (ps) 

Fig.  11  -  Time-History  From  Wall-Mounted  Transducers  for  PDC-M30  (Problem  #4);  210  m/s  Impact 
on  60.5%  TMD  VVC-I40.  Comparison  of  Workshop  Predictions  From  Weston  with  Data  from  Sandusky 
et.  al.  (3|.  Gage  #1  is  Located  at  38.1  mm,  and  Gfige  #2  at  76.2  mm;  Tube  Length  is  146.8  mm. 
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possibility  that  the  real  compactioa  wave  is  reactive  and  exhibits  rate-dependent 
material  behavior. 

Since  a  number  of  the  modelers  are  assuming  that  the  compaction  wave  is 
not  inert,  there  were  additional  discussions  about  what  reaction  is  being  modeled. 
There  appear  to  be  two  broad  categories  here,  (a)  condensed-phase  "hot  spot" 
reactions,  and  (b)  gas-phase  reactions.  Some  models  include  both.  There  was  no 
consensus  on  how  to  model  hot  spot  reactions.  A  couple  of  the  models  which 
include  gas-phase  reactions  assume  that  the  reactants  originate  as  products  from 
the  condensed-phase  reactions.  Postulating  the  existence  of  reactive  gas-phase 
species  leads  to  an  interesting  question  about  the  extent  of  heat  release  in  the 
condensed-phase  reactions.  The  early  PDC  experiments  on  TS-3659  in  a  Lexan 
tube  showed  no  indication  of  visible  radiation,  which  would  seem  to  constrain 
gas-phase  temperatures  well  below  an  equilibrium  flame  temperature.  Thus  the 
sketchy  evidence  so  far  would  suggest  the  initial  gas-phase  reactants,  if  they 
exist,  are  intermediate  products  which  then  could  react  further  at  some  later 
time.  Release  of  this  final  heat  of  combustion  would  be  a  contributor  to  ignition 
of  the  remaining  portion  of  the  compressed  aggregate. 

The  results  of  the  Workshop  might  be  summarized  with  the  following; 

(1)  The  compaction  front  locus  (position  vs.  time)  and  hence  compaction  wave 
speed  are  forgiving  parameters.  Agreement  with  this  data  does  not  guarantee 
that  the  model  can  predict  the  pressure  transducer  records. 

(2)  The  impact-generated  compaction  wave  apparently  sets  off  an  intense 
competition  in  the  compressed  aggregate  between  various  reaction  sources  and  heat 
loss;  the  outcome  determines  the  time  of  a  runaway  pressure  field.  Current 
models  have  difficulty  making  this  prediction.  Besides  the  numerical  solution 
complexities,  the  primary  reason  is  uno'rtainty  about  the  reaction  sources,  how  to 
model  these  sources  and  then  calibrate  the  inevitable  unknown  parameters. 

(3)  When  comparing  model  predictions  with  the  PDC  experimental  data: 

(a)  several  of  the  models  were  able  to  predict  the  general  wave  mechanics 
(time  at  which  each  gage  reports,  wave  thickness,  etc.) 

(b)  the  wall-mounted  transducer  data  show  an  almost  linear  "decay"  after 
the  initial  report  value,  until  the  onset  of  significant  reaction.  None  of  the 
models  predict  this,  and  it  has  not  been  explained. 

(c)  the  initial  value  reported  by  the  second  transducer  is  less  than  that 
reported  by  the  first  transducer.  Wall-boundary  friction  is  the  likely  cause. 
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APPENDIX  -  SCHEMATIC  DESCRIPTION  OF  MODELS 


33 


BAER  (SNLA)  [see  Ref.  13] 

0  Compressible  solid  phase, 
o  Zero  wall  friction. 

o  Equilibrium  stress  state,  0  ,  is  fit  to  the  combined  dynamic  data  (PDC)  and 
quasi-static  compaction  c^ta  (Elban,  et.  al.) 

o  Does  not  account  for  stress  deviators  (each  phase  represented  with  a  pressure). 

o  Finite-rate  porosity  change  based  on  "compaction  viscosity"  -  leads  to 
prediction  of  finite-thickness  compaction  wave. 


Mass  Source  Term 


Compressive  Reaction 


O 

T  X 


H 


+  Grain  Burning 

^  0  only  when  T.  >  T 
heat  transfer  to  cold  solid 
determines  Tj(t) 


Compressive  reaction  - 

2 

o  1  /  Tr  =  constant  * 

o  Energy  released  in  reaction  depends  upon  delay  time  "  I  ". 

AE  =  Q*  +  (  AEp  -  Q*)  H(I-l)  where  H(I-l)  =  1  when  I  >  1. 
0  Delay  time  "  I  "  is  tracked  along  gas-phase  streamlines,  such  that 
D(  I  )/D(tg)  =  1  /  where  =  constant  *  P^^ 

Pressure  or  Stress 
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HSIEH  &  i™  (NSWC)  [see  Ref.  14) 
o  Incompressible  solid  phase, 
o  Zero  wall  friction. 

o  Equilibrium  stress  state  is  fit  to  qu2isi-static  compaction  data. 

o  Does  not  account  for  stress  deviators  (hydrostatic  solid  stress  and  gas  pressure 
only). 

o  No  rate-dependent  resistance  to  porosity  change  ("compaction  viscosity"  =  0). 

o  Drag  force  between  solid  and  gas  phases  was  neglected  in  Workshop 
computations. 

Uses  GOUGH'S  "XNOVAK"  Code - with  modifications 


Mass  Source  Term  =  Grain  Burning 

o  Single  step  reaction  releases  full  chemical  energy  of  solid  phase 
o  Onset  of  this  reaction  is  delayed  until  achieving  f  ,  where 

*  2 

f  =  1  r.  dt,  for  each  Eulerian  control  volume 

Jo 

*  2  * 
the  value  of  f  follows  from  Sandusky's  correlation,  Tj  At  =  constant  =  f  . 


Pressure  or  Stress 
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KOOKER  (BRL)  [see  Ref.  15] 

0  Compressible  solid  phase. 

o  Zero  wall  friction. 

o  Equilibrium  stress  state  is  fit  to  quasi-static  compaction  data. 

o  No  explicit  use  of  stress  de viators  (each  phase  represented  with  a  pressure). 

o  Instantaneous  porosity  adjustment  (compaction  viscosity  =  0)  to  equilibrium 
stress  state. 

o  Compaction  waves  are  shock  waves  with  zero  thickness  (method-of- 
characteristics  solution) 

o  Compaction  wave  front  induces  small  amount  of  reaction  (produces  gas-phase 
intermediate  products) 


Mass  Source  Term  =  Gas-Phase  Reaction  +  Grain  Burning 

^  0  only  when  T.  >  T 
heat  transfer  to  cold  solid 
determines  Tj(t) 

Gas-phase  reaction  "begins"  with  intermediate  products  from  wave-induced  reaction. 

o  Reaction  is  delayed  until  tj^  <  0 

D(t.J/D(t„)  =  -1,  where  t.  (0)  =  t.  °  comes  from  P  ^At  =  constant 
ig"  g'  ig'  ’  ig  m 

o  Amount  of  energy  released  is  a  function  of  gas-phase  temperature, 

D(Ae  )/D(t  )  =  constant  *  exp(-constant  /  T  ) 

s  s 

Pressure  or  Stress 
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Distance 


PRICE  (NWC)  [see  Ref.  16] 


0  Compressible  solid  phase. 

0  Accounts  for  wall  friction. 

0  Equilibrium  stress  state  is  fit  to  quasi-static  compaction  data. 

0  Mixture  stress  tensor  is  a  pressure,  but  intragranular  stress  is  function  of 

parameters  called  shear  modulus,  bulk  modulus,  etc.  [At  Workshop,  reported 
radial  stress  =  0.6  *  mixture  pressure.] 

0  No  rate-dependent  resistance  to  porosity  change  ("compaction  viscosity"  =  0) 
other  than  inertial  forces. 

Heavily  Modified  Version  of  Joint  Venture  Two-Phase  Flow  Code  ( —  roots  to 

Krier,  et.  al.) 

No  Artificial  Distinction  Between  Ignition  Transient  Combustion 
All  Reaction  &;  Energy  Release  are  Incorporated  into  "Two-Path  Kinetics" 
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WESTON  (Brobeck/LLNL)  [see  Ref.  17] 

0  Compressible  solid  phase, 

o  Accounts  for  wall  friction. 

o  Equilibrium  stress  state  is  fit  to  quasi-static  compaction  data, 
o  Extensive  analysis  based  on  stress  deviators. 

o  No  rate-dependent  resistance  to  porosity  change  ("compaction  viscosity"  =  0). 


Modified  "R-DUCT"  Model 


Mass  Source  Term  ==  Gas-Phase  Reaction 


Gas-Phase  Reaction  - 


+ 


Grain  Burning 


0  only  when  T.  >  T 
heat  transfer  to  coid  solid 
determines  Tj(t) 


o  Reaction  releases  full  chemical  energy  of  solid  phase. 

o  Delay  time  before  onset  of  reaction  is  controlled  by  parameters  fit  to 
"ODTX"  data. 

plus 

’  0  An  additional  reaction  in  the  leading  compaction  front  generates  some  gas, 
according  to 

Production  rate  =  constant  *  strength  of  wave  *  compaction  rate  at 

local  cell 


Pressure  or  Stress 
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ii.nciass.,  unlimited)  12  Administrator 

it  nciaa,  limited)  2  Dcfensc  Tcchiucal  Info  Center 

(Ciassincd)  2  ATTN:  DTIC-DDA 

Cameron  Station 
Alexandria,  VA  22304-6145 


1  HQDA  (SARD-TR) 

WASH,  DC  20310-0001 

1  Commander 

US  Army  Materiel  Command 
ATTN;  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 


1  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOC) 

Redstone  Arsenal,  AL  35898-5010 

1  Commander 

US  Army  Tank  Automotive  Command 
ATTN:  AMSTA-TSL  (Technical  Liteary) 
Warren,  MI  48397-5000 

1  Director 

US  Army  TRADOC  Analysis  Command 
ATTN:  ATAA-SL 

White  Sands  Missile  Range,  NM  88002-5502 


1  Commander 

US  Army  Laboratory  Command 
ATTN;  AMSLC-DL 
Adelphi,  MD  20783-1145 


(ciiu.  only)  1  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 
Fort  Benning,  GA  31905-5660 


2  Commander 

Armament  RD&E  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-MSI 
Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

Armament  RD&E  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-TDC 
Picatinny  Arsenal,  NJ  07806-5(X)0 

1  Director 

Benet  Weapons  Laboratory 
Armament  RD&E  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-LCB-TL 
Watervliet,  NY  12189-4050 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCARESP-L 
Rock  Island,  IL  61299-5000 


(UbcIim.  only)  1  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 

1  The  Rand  Corporation 
P.O.  Box  2138 

Santa  Monica,  CA  90401-2138 

(Clan,  mly) 

1  Air  Force  Armament  Laboratory 
4TTN-  AFATf/nLODL 
Eglin  AFB,  FL  32542-5000 

Aberdeen  Provina  Ground 
Dir.  USAMSAA 

ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 

Cdr,  USATECOM 

ATTN:  AMSTE-TO-F 
Cdr,  CRDEC,  AMCCOM 

ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSl 


1  Commander 

US  Army  Aviation  Systems  Command 
ATTN:  AMSAV-DACL 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120-1798 


1  Director 

US  Army  Aviation  Research 
and  Technology  Activity 
Ames  Research  Center 
Moffeu  Field,  CA  94035-1099 
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Copies  Or^anlzat Ion 

1  Commander 

Armament  R&D  Center 
ATTN:  SMCAR-SCA-T, 

L.  Stiefel 
Bids.  455 

Dover,  NJ  07801-0001 

1  Commander 

ns  R&D  Center 
ATTN:  SMCAR-LCA-G, 

A  ..I .  Beardell 
Dover,  NJ  07801-0001 

1  Commander 

Armament  R&D  Center 
Scientific  &  Tech  Div 
ATTN:  SMCAR-TSS 

Rldg.  59 

Dover,  NJ  07801-0001 

1  Director 

Army  Def  Ammo  Ctr  &  Sch 
ATTN:  SMCAC-ASA  (IIRC)  , 

A.H.  Bausman 
Savanna,  IL  61074 

1  Commander 

Army  Fgn  Sci  &  Tech  Ctr 
ATTN:  CH,  H.M.  Sheinfeld 

220  7th  Street,  NE 
Charlottesville,  VA  22901 

1  Commander 

Army  Laboratory  Command 
Mat  Tech  T^b 

ATTN:  SLCMT-TML-Tech  Lib, 

M. M.  Murphy 
Watertovm,  MA  02172 

1  Commander 

Army  Mat  Command 
ATTN:  AMCSF,  Safety  Office 
5001  Elsenhower  Avenue 
Alexandria,  VA  22333 


No.  Of 

Copies  Organization 

1  Commander 

Army  Missile  Command 
Redstone  Sci  Inf  Ctr 
ATTN:  DOC  SEC 
Bldg.  4484 

Redstone  Arsenal,  AL  35898-5000 

2  Commander 

Army  Msl  Command 
Redstone  Sci  Inf  Ctr 
ATTN:  AMSMl-RD-PR-S, 

S.P.D.  Smith 
AMSMI-XO,  F.F. .  Hart 
Redstone  Arsenal,  AL  35898-5000 

2  Director 

(IS  Army  Research  Office 
ATTN:  DRXRO-PP-LIB 

DEXTR-PP-LIR,  D.M.  Mann 
P.O.  Box  12211 
Research  Triangle  Park,  NC 
27709 

2  Commander 

Army  Strategic  Defense 
Command  -  Huntsville 
ATTN:  DASI>H-MPL,  A.B.  Dumas 

BMDSC-AD-LIB, 

D.C.  Sayles 
P.O.  Box  1500 
Huntsville,  AL  35807 

1  Commander 

White  Sands  Missile  Rnage 

ATTN:  Tech  Lib 

l^ite  Sands  Missile  Range, 

NM  88002 

3  Commander 

Naval  Air  Sys  Command 
ATTN:  NAIR-32/931A, 

C.E.  Kovalenko 
ATR  5004, 

S.N.  Block  (54042) 
AIR-5004- Tech  Lib 
Washlngti)!!,  DC  20361“9300 
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Copies  Organization  Copies  Organlzat Ion 


1  Commanding  Officer 

Naval  Explo  Ord  Dspl  Fac 
ATTN;  Lib  Dlv 

Indian  Head,  MD  20640 

2  Commanding  Officer 

Naval  Tech  Intelligence  Ctr 

ATTN:  Code  DS-311 

4301  Sultland  Road 
Washington,  DC  20395-5020 

3  Commanding  Officer 

Naval  Ordanance  Station 
ATTN:  Code  5253,  T.C.  Smith 

Code  5246,  Tech  Lib 
Code  OE-5 ,  J.A.  Hannum 
Indian  Head,  MD  20640 

1  Superintendent 

Naval  Postgraduate  School 
ATTN:  Code  1424,  Libs  Dir 

Monterey,  CA  93943 

1  Chief 

Naval  Research 
ATTN:  Code  1132P, 

R.S.  Miller 
800  N.  Quincy  Street 
Arlington,  VA  22217 

1  Commanding  Officer 

Naval  Research  Lab 
ATTN;  Code  2627 
Washington,  DC  20375 

1  Commander 

Naval  Sea  Systems  Command 
ATTN;  SEA-06 H, 

M.R.  Vanslyke 
Washington,  DC  20362 

1  Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA-09 B312,  Tech  Lib 

National  Center  Bldg.  3 
Washington,  DC  20362 


2  Commanding  Officer 

Naval  Ship  Weapon  System 

Engineering  Station 
ATTN:  Code  5125, 

E.C.  Kalsteln 
Code  4V21 
M.W.  Walker 

Port  Hueneme,  CA  93043 

2  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  G-20,  J.L.  East 

Code  E231,  Tech  Lib 
Dahlgren,  VA  22448 

2  Commander 

Naval  Surface  Weapons  Centei 
White  Oak  Laboratory 
ATTN:  Code  K-22 , 

R.J.  Edwards 
Code  E432,  S.  Happel 
10901  New  Hampshire  Avenue 
Silver  Spring,  MD  20903 

7  Commander 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATTN:  Code  R-13 

R. R.  Bernecker 

S.  Jacobs 
H.W.  Sandusky 
D.  Price 

B.  Clancy 
K.  Kira 

C. W.  Dickinson 
10901  New  Hampshire  Avenue 
Silver  Spring,  MD  20903 

1  Commanding  Officer 

Naval  Underwater  Sys  Cmd 
ATTN;  021312,  Tech  Lib 
Newport,  RT  02840 
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C  o  o i e  s  Org anlzatlon 

4  Cornnander 

NavaL  Weapons  Center 
■\TTrJ;  Code  MTIH, 

A .0 .  Smith 

Code  3272,  R.W.  Feist 
Code  327B,  R.F.  Vetter 
Code  3431 

China  Lake,  CA  93553 

4  Commander 

Naval  Weapons  Center 
ATTN:  Code  3891, 

A.  Atwood 
H.P.  Richter 
C.  Price 
F.  Lundstrom 
China  Lake,  CA  93555 

1  Commanding  Officer 

Naval  Weapons  Station 
ATTN:  Code  322, 

T.W.  Lohrmann 
Concord,  CA  94520 

2  Director 

Navy  Strat  Sys  Proj  Ofc 

ATTN:  SP-2731,  J.E.  Culver 

Tech  Lib  Br  Hd 
Washington,  DC  20376 

1  Commander 

Operations  T&E  Force  Atlantic 
ATTN:  Tech  Lib 

Norfolk,  VA  23511 

1  Commander 

AAMRL 
ATTN:  TH 

Wright-Patterson  AFB, 

OH  45433 

1  Commander 

AD 

ATTN:  ENMP,  R.H.  Rhode 
Eglin  AFB,  FL  32542 


No.  Of 

Copies 

Organization 

1 

Commander 

AEDC 

ATTN:  Tech  Lib 

Arnold  AFS,  TN 

37389 

1 

Commander 

AFFTC 

ATTN:  Tech  Lib 

Stop  238 

Edwards  AFB,  CA 

93523 

10 

Commander 

AFRPL 

ATTN:  DYC,  D.P. 

Weaver 

MKPL,  C.E.  Merrill 
MKPB,  D.T.  Thrasher 
MKPA,  J.L.  Trout 
MKAS,  J.H.  Clark 
LKCP,  J.W.  Marshall 
DY,  R.L.  Geisler 
MKA 
MKP 

Tech  Lib 

Edwards  AFB,  CA  93523 

1  Command!  .. 

AFSC 

ATTN:  DLX-T,  R.E.  Smith 

Andrews  AFB 
Washington,  DC  20334 

1  Commander 

AFWAL  (GLISL) 

ATTN:  Tech  Rpts  Sec 

Wright-Patterson  AFB, 

OH  45433 


1  Commander 

ASD 

ATTN:  YYEF,  .1  .v  .  Braum 

Wright-Patterson  AFB, 

OH  45433 
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No.  Of 
Copies 

Organization 

Commander 

1 

Purdue  University 

ESMC 

Aeronautics  &  Astronautics 

ATTN:  ETR/RA,  M.  Brown 

ATTN:  J.R.  Osborne 

Patrick  AFB,  FL  32925 

Grissom  Hall 

Commander 

FTD 

1 

West  Lafayette,  IN  47907 

Southwest  Research  Institute 

ATTN:  TQTA,  A.  Crowder 

ATTN:  H.J.  Cry ting 

Wright-Patterson,  AFB 

P.O.  Drawer  28510 

OH  45433 

San  Antonio,  TX  78284 

Commander 

FTD 

ATTN:  SDBP 

Wright-Patterson  AFB, 

OH  45433 

Commander 
Ogden  Ale 

ATTN:  jMANPA,  A.J.  Inverse 

Bldg.  1941 

Hill  AFB,  UT  84056 

Commander 

SAC 

ATTN:  NRI-STINFO 
Offutt  AFB,  NB  68113 

California  Inst  Technology 
Jet  Propulsion  Lab,  Lib  Op  Gp 
ATTN:  Lib  Acqs/ Standing 

Orders 

4800  Oak  Grove  Drive 
Pasadena,  CA  91103 

IIT  Rsch  Inst 
ATTN:  Doc  Lib 

10  W.  35th  Street 
Chicago,  IL  60616 

Pennsylvania  State  Unlv. 
Applied  Rsch  Lab 
ATTN;  K.K.  Kuo 
312  Mechanical  Eng.  Bldg. 

P.O.  Box  30 

State  College,  PA  16801 


4  Lawrence  Livermore  National 

Laboratory 
ATTN:  L.  Lee,  L324 

M.  Costantino,  L324 

A.  Weston,  L324 
E.  James,  L324 

P.O.  Box  808 
Livermore,  CA  94550 

6  University  of  California 

Los  Alamos  National  Lab 
ATTN;  Report  Lib,  R.  Rabie 
J.  McAfee,  J960 

B.  Asay,  J960 
J.  Dick,  P952 
J.  Ramsay,  J960 

C.  Forest 
P.O.  Box  1663 

Los  Alamos,  NM  87545 

1  University  of  Illinois  at 

U^b?;^a 

ATTN:  H.  Krier 
144  ME  Building 
Urbana,  IL  61801 

1  University  of  Illinois  at 

Urbana 

Engineering  Exp  Station 
ATTN:  Doc  Lib 

1308  W.  Green  Street 
Urbana.  IL  61801 

1  AAI  Corporation 

ATTN:  Tech  Lib 

P.O.  Box  126 
Hunt  Valley,  MD  21030 


No .  Of 

Copies  Organization 

1  Aerojet  General  Corp 

Sacramento  Facility 
ATTN:  Tech  Inf  Ctr 

P.O.  Box  13222 
Sacramento,  CA  95813 

1  Aerojet  Tactical  Systems 

ATTN:  R.  Mlronenko 

P.O.  Box  13400 
Sacramento,  CA  95813 

1  Aerospace  Corp 

ATTN:  S.B.  Crowe 

P.O.  Box  92957 

Los  Angeles,  CA  90009 

1  Air  Products  &  Chem,  Inc. 

ATTN:  C.R.  Anderson 

P.O.  Box  538 
Allentown,  PA  18105 

1  ARINC  Research  Corp 

Annapolis  Sci  Ctr 
ATTN:  R.  Blair 
2551  Riva  Road 
Annapolis,  MD  21401 

1  Atlantic  Research  Corp 

ATTN:  Tech  Inf  Ctr 

7511  Wellington  Road 
Gainesville,  VA  22065 

1  Boeing  Company 

Kent  Library 

ATTN:  8K-38  Lib  Proc  Supv 
Mail  Stop  74-60 
P.O.  Box  3707 
Seattle,  WA  98124 

1  Cal span  Corp 

ATTN:  V.M.  Young 
P  .0.  Box  400 
Buffalo,  NY  14225 

1  Ford  Aerospace  and  Comra 

Aeronutronic  Div 
ATTN:  Tech  Inf  Svc,  DDC  ACQS 
Ford  &  Jamboree  Roads 
Newport  Beach,  CA  92658 


No.  Of 

Copies  Organization 

1  General  Dynamics  Corp 

Pomona  Division 

ATTN:  Div  Lib,  MZ-4-20-CH  LIB 
Mail  Zone  6-20 
P  .0  .  Box  2507 
Pomona,  CA  91766 

1  General  Dynamics  Corp 

Convair  Div 
ATTN;  U.J.  Sweeney 
P.O.  Box  85357 
San  Diego,  CA  92138 

1  General  Elec  Co 

Armament  &  Elec  Sys  Dept 
ATTN:  Eng  Lib 
Lakeside  Avenue 
Burlington,  VT  05402 

1  Grumman  Aerospace  Corp 

Engineering  Lab 
ATTN:  H.B.  Smith 
MS  iJOl-35 

S.  Oyster  Bay  Road 
Bethpage,  LI,  NY  11714 

1  Hercules,  Inc. 

Aerospace  Division 
Allegany  Ballistics  Lab 
ATTN:  LIB 

P.O.  Box  210 
Cumberland,  MD  21502 

2  Hercules,  Inc. 

Bacchus  Works 

ATTN:  Library  M/ S  H, 

C.F.  Partridge 
MS-J12D, 

G.  Butcher 
P.O.  Box  98 

Magna,  UT  84044-0098 

1  Hercules,  Inc. 

ATTN;  Tech  Inf  Ctr, 

Bldg.  A- 100 
P.O.  Box  548 
McGregor,  TX  76657 
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1  HI,  Shear,  Corp 

ACHTech  Div 

.ATTN:  Eng  Dir,  D.E.  Olander 

2600  Skypark  Drive 
Torrance,  CA  90509 

1  Hughes  Aircraft  Co 

Electro  Optical  &  Data  Sys  Gp 

ATTN:  B.W.  Campbell 

Bldg.  El  El  10 

P.O.  Box  902 

El  Segundo,  CA  90245 

3  Lockheed  Msl  &  Space  Co 

Tech  Inf  Ctr 

ATTN:  Reports  Org  52-52,  B.  201, 
Lockheed  Huntsville 
Facility 

ORG  83-10,  Bldg.  157-3W, 

C.  Havlik 

Advanced  Aeronautics, 

T.  Harsha 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

1  Ltv  Aerospace  &  Def  Co 

ATTN:  LIB  3-58200 
Mail  Stop  EM-08 
P.O.  Box  655907 
Dallas,  TX  75265-5907 

1  Marquardt  Co 

Sub  Tsc  Def  &  Space  Gp  Inc 

ATTN :  LIB 

P.O.  Box  9104 

Van  Nuys,  CA  91409 

1  Martin  Marietta  Corp 

Orlando  Aerospace 
ATTN:  MP-30  -  TIC 
P.O.  Box  5837 
Orlando,  FL  32855 

1  Martin  Marietta  Corp 

Denver  Aerospace 
ATTN:  Rsch  Lib  2655 
P.O.  Box  179 
Denver,  CO  80201 


1  McDonnell  Douglas  Corp 

ATTN:  R.D.  Detrick 

P.O.  Box  516 
St.  Louis,  MO  63166 

1  Morton  Thiokol,  Inc. 

Tech  Staff,  Aerospace  Group 
ATTN:  T.F.  Davidson 

3340  Airport  Road 
Ogden,  UT  84405 

1  Morton  Thiokol,  Inc. 

Elkton  Div,  Aerospace  Group 
ATTN:  D.J.  McDaniel 

P.O.  Box  241 
Elkton,  MD  21921 

1  Morton  Thiokol,  Inc. 

Wasatch  Div,  Aerospace  Group 
ATTN:  J.E.  Hansen 
P.O.  Box  524 

Brigham  City,  UT  84302 

1  Morton  Thiokol,  Inc. 

Huntsville  Div 
ATTN;  H.H.  Sellers 
Huntsville,  AL  35807 

1  Morton  Thiokol,  Inc. 

Longhorn  Div 
ATTN:  Tech  Lib 

P.O.  Box  1149 
Marshall,  TX  75670 

1  Olin  Corp 

Ordnance  Prods 
ATTN:  M.R.  Rhine 

P  .0.  Box  278 
Marion,  IL  62959 

1  Raytheon  Co 

Missile  Sys  Div 
ATTN;  J.  Durben 
Hartwell  Road 
Bedford,  MA  01730 
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1  Rockwell  International 

Rocketdyne  Div 
ATTN:  Tech  Inf  Ctr 

6633  Canoga  Avenue 
Canoga  Park,  Ca  91304 

I  Rockwell  International 

Space  Sys  Gp 

ATTN:  TIC-D201-300-AJ01 

12214  Lakewood  Blvd . 

Downey,  CA  90241 

3  Sandia  National  Labs 

Mail  Svc  Sec 
ATTN:  Lib  3141, 

J .  Van  Berkel 
M.R.  Baer,  MS  1513 
J.  Cummings,  MS  1513 
P.O.  Box  5800 
Albuquerque,  NM  87185 

1  Stone  Engineering  Company 

ATTN:  Scty  Off,  D.  Owens 
805  Madison  Street  ( STE  2C) 
Huntsville,  AL  35801 

1  Sundstrand  Corp 

ATTN:  C.A.  Hull 
P.O.  Box  7003 
4751  Harrison  Avenue 
Rockford,  IL  61125 

1  Talley  Inds  Arizona  Inc 

ATTN:  K.  St.  Clair 

P.O.  Box  849 
Mesa,  AZ  85201 

1  Textron,  Inc. 

AVCO  Rsch  Lab  Textron 
ATTN:  L.T.  Nazzaro 
2385  Revere  Beach  Pkwy 
Everett,  MA  02149 

1  Textron  Inc 

Bell  Aerospace  Co  Div 
ATTN:  L.M.  Breslauer 

P.O.  Box  1 
Buffalo,  NY  14240 


5  TRW  Inc 

TRW  Electronics  &  Defense 
Sector 

ATTN:  R.C.  Reeve,  San 
Bernardino 
R.J.  Christiansen, 
01/2010 
S/1930 
C.M.  Orr 

R.  Fleeter,  Rl/1028 
One  Space  Park 
Redondo  Beach,  CA  90278 

1  United  Technologies  Corp 

Chemical  Systems  Div 
ATTN:  Tech  Lib 

P.O.  Box  50015 
San  Jose,  CA  95150-0015 

1  United  Technologies  Corp 

Research  Center 
ATTN:  M.E.  Donnelly 
400  Main  Street 
East  Hartford,  CT  06108 

1  Universal  Propulsion  Co 
ATTN:  H.J.  McSpadden 
Black  Canyon  Stage  1 
Box  1140 

Phoenix,  AZ  85029 

2  Natl  Aeronautics  &  Space 

Admin 

George  C.  Marshall  Space 
Flight  Center 
ATTN:  SA42,  J.Q.  Miller 
AS24D 

Marshall  Space  Flight  Center, 
AL  35812 

3  Natl  Aeronautics  &  Space 

Admin 

Office  of  Management 
ATTN:  Admin  Div  NHB-12, 

Lib  Sec 

RTP-6 ,  F.W.  Stephenson 
Code  DS,  W.R.  Frazier 
Washington,  DC  20546 
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1  Natl  Aeronautics  &  Space 

Admin 

Wallops  Flight  Center 

ATTN:  Code  250.9,  J.N.  Foster 

Wallops  Island,  WA  23337 

1  Natl  Aeronautics  &  Space 

Admin 

Langley  Research  Center 
ATTN:  MS-185,  Tech  Lib 
Hampton,  VA  23665 

1  Natl  Aeronautics  &  Space 

Admin 

Lewis  Research  Center 
ATTN:  D.  Morris 

21000  Brookpark  Road 
Cleveland,  OH  44135 

1  Natl  Aeronautics  &  Space 

Admin 

Scientific  Tech  Inf  Fac 
ATTN:  Accessioning  Dept 
P.O.  Box  8757 
Baltlmore-Washington  Inti 
Airport,  MD  21240 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
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